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Abstract 

Wide-band Suzaku data on the merging cluster Abell 3667 were examined for hard X-ray emission in 
excess to the known thermal component. Suzaku detected X-ray signals in the wide energy band from 
0.5 to 40 keV. The hard X-ray (> 10 keV) flux observed by the HXD around the cluster center cannot 
be explained by a simple extension of the thermal emission with average temperature of ~ 7 keV. The 
emission is most likely an emission from a very hot (kT > 13.2 keV) thermal component around the cluster 
center, produced via a strong heating process in the merger. In the north-west radio relic, no signature of 
non-thermal emission was observed. Using the HXD. the overall upper-limit flux within a 34' x 34' field-of- 
view around the relic is derived to be 5.3x 10~ 12 erg s _1 cm -2 in the 10-40 keV band, after subtracting the 
ICM contribution estimated using the XIS or the XMM-Newton spectra. Directly on the relic region, the 
upper limit is further tightened by the XIS data to be less than 7.3 x 10~ 13 erg s _1 cm~ 2 , when converted 
into the 10-40 keV band. The latter value suggest that the average magnetic field within the relic is higher 
than 1.6 fiG. The non-thermal pressure due to magnetic fields and relativistic electrons may be as large as 
~ 20% of the thermal pressure in the region. 

Key words: X-rays: galaxies: clusters: individual (Abell 3667) — galaxies: magnetic fields - 
radiation mechanisms: non-thermal 



1. Introduction 

Enigmatic extended radio sources in rich galaxy clus- 
ters have been known for over 30 years (Willson 1970; for 
recent review, see Feretti 2005). Sources located at the 
cluster center are referred to as "radio halos" , while those 
on the cluster periphery are called "relics" . Steep spectra 
and spectral cutoffs at a few GHz detected in several ha- 
los/relics indicate that relativistic electrons with typically 



GeV energy are undergoing cooling due to synchrotron 
and inverse Compton (IC) emission. In every case, the 
radio halos / relics are found in irregular clusters which are 
apparently undergoing mergers. This suggests that the 
radio emitting electrons arc accelerated by shocks or tur- 
bulence generated by energetic cluster merging. 

The same relativistic electrons scatter Cosmic 
Microwave Background (CMB) photons up to the 
hard X-ray band (IC emission). By comparing the radio 
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and the hard X-ray emissions, not only the electron 
energy density but also the intra-cluster magnetic field 
intensity can be estimated (e.g. Sarazin 1988). Detecting 
the IC emission is difficult because of the strong ICM 
component dominant below ~ 20 keV and the low 
expected flux of the emission in hard X-rays above this 
energy. Currently, there are only a few reports claiming 
its detection, mostly from the Beppo-SAX PDS in the 
energy band around 40 keV or higher (e.g. Fusco-Femiano 
et al. 1999; Nevalainen et al. 2004), though a few of the 
detections are still controversial (Rossetti & Molendi 
2004; Fusco-Femiano et al. 2007). There are reports from 
RXTE (e.g. Petrosian et al. 2006; Rephaeli et al. 2006) 
as well. Another case includes the ASCA GIS results 
around 4 keV on galaxy groups (Fukazawa et al. 2001; 
Nakazawa et al. 2006). These results would require a 
magnetic field as low as 0.1 /zG (e.g. Fusco-Femiano et 
al. 1999), which at the first glance is inconsistent with 
the radio rotation measure estimation of 2 ~ 10 /xG 
(e.g. Clarke ct al. 2001; Carilli & Taylor 2002). Further 
observations by independent instruments of IC and other 
hard X-ray emission processes (e.g. Inoue et al. 2005), is 
thus desirable. 

Significant heating of the ICM itself must also be taking 
place in merging clusters (e.g. Takizawa 2000; Ricker & 
Sarazin 2001). X-ray hardness ratio images sometimes 
show temperatures exceeding 10 keV (e.g. Watanabe et 
al. 1999; Briel et al. 2004). Since the emission can be 
multi-phase within the line of sight, the actual highest 
temperature can be much higher, though it is not easy to 
identify using the contemporary imaging X-ray detectors 
working below < 12 keV. 

In many ways, Abell 3667 is the ideal cluster to study 
mergers, radio relics and hard X-ray emission. It is a very 
bright X-ray cluster at a low redshift of z = 0.0556 (Strublc 
& Rood 1999). The ROSAT and ASCA observations 
showed that it is a spectacular merger with shock heated 
gas (Markevitch et al. 1999). Chandra and XMM-Newton 
observations have shown much evidence for an ongoing 
merger, such as a cold front (Vikhlinin et al. 2001a, b) and 
highly inhomogeneous temperature structure with tem- 
perature ranging from 4 to > 10 keV (Mazzotta et al. 2002; 
Briel et al. 2004). The optical galaxy distribution shows 
elongation towards the north-west south-east axis (Sodrc 
et al. 1992; Johnston-Hollitt et al. 2008), supporting the 
binary merger scenario as well. 

The cluster is famous for its pair of radio relics (e.g. 
Roettgering et al. 1997). The relic to the north-west is 
the brightest and largest among the diffuse radio sources 
associated with cluster of galaxies, having a flux of 3.7 
Jy ( Johnston- Hollit 2004) and a width of ~ 20' at 1.4 
GHz. It is located about 30' or 2 Mpc from the clus- 
ter X-ray centroid. The radio signal is detected from 85 
MHz to 2.3 GHz, and its average photon index is T = 2.1 
(Roettgering et al. 1997). Here the photon index T is 
defined as A^ p h ton(£') = No x E~ T . The radio relics have 
very sharp outer edges and the radio spectra steepen with 
distance from the edge. Thus, the radio relics are consid- 
ered to reflect the position of the particle acceleration in 



a merger shock. The cluster is also a good candidate for 
detection of IC emission. The Beppo-SAX PDS provided 
a rather marginal evidence of hard excess emission from 
the cluster (Fusco-Femiano et al. 2001; Nevalainen ct al. 
2004). However, the large field of view (FOV) of the PDS 
makes the data strongly affected by contamination from 
thermal emission from the entire cluster, mainly from its 
central region, and possible AGNs. 

The hard X-ray detector (HXD; Takahashi et al. 2006; 
Kokubun et al. 2006) onboard Suzaku (Mitsuda et al. 
2006) is characterized by its low detector background in 
the 10-40 keV band and its narrow FOV of 34' full width 
at half maximum (FWHM). Using the HXD, wc can spa- 
tially distinguish hard X-ray components from the cluster 
center and the north-west relic, with the highest sensitiv- 
ity in the 10-40 keV band. In addition, the X-ray CCD 
cameras (XISs; Koyama et al. 2006) onboard Suzaku are 
characterized by their large effective areas and low and 
stable detector backgrounds, which make them very pow- 
erful devices to observe extended hard emission. 

In May 2006, we observed Abell 3667 cluster using 
Suzaku with three pointings running from the north-west 
relic to the center, each separated by ~ 17'. Here we report 
the results obtained by these observations. Observation 
logs and data reduction are described in the next section. 
Section 3 is devoted to the spectral and imaging analysis, 
followed by the discussion in section 4 and conclusions in 
section 5. In this paper, we utilized canonical cosmology 
value of H = 70 Mpc km^s. At the redshift of z = 0.0556, 
1' corresponds to 62 kpc. For a mean cluster temperature 
of 7.2 keV, the virial radius should be r 2 oo ~ 2.6 Mpc = 42' 
(e.g. Neumann & Arnaud 1999). Unless otherwise noted, 
all errors are at the 90% confidence limit. 

2. Observation and Data Reduction 

2.1. Observation logs 



Table 1. Log of the three Abell 3667 observations. 



Obs. 


Start (UT) 


End (UT) 


Exp. (ks) * 


NW.Relic 
NW.17off 
Center 


05/03 17:47 
05/06 07:03 
05/06 17:40 


05/06 07:02 
05/06 17:39 
05/07 09:13 


81.2/45.6 
16.7/11.4 
20.5/7.2 



* The XIS/HXD exposure time after the screening as presented in 
the text. 



The three Abell 3667 pointing observations of Suzaku 
were carried out successively in the beginnng of May 2006. 
As already noted, individual pointing positions are offset 
by ~ 17', so that the first observation covers the north- 
west relic, and the last one on the cluster center, with 
continuous XIS coverage. Hereafter, we call these obser- 
vations, the NWR, the 17' offset and the center pointing, 
respectively. Details of the observation time are summa- 
rized in table 1. The FOVs of individual pointings com- 
pared with the X-ray sky image of PSPC are also pre- 
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sented in appendix 3. 

The XIS was operated in the normal timing and full 
window mode. All the 64 PINs of the HXD were operated 
in the nominal bias voltage of 500 V. All the data was 
processed with the Suzaku pipeline processing of version 
2.0.6.13. 

2.2. Data reduction of the XIS 

2.2.1. Screening and background estimation 

The XIS data were processed via default screening cri- 
teria. Events with a GRADE of 0, 2, 3, 4, 6 and STATUS 
< 1024 were extracted. Data screening criteria with HK 
files are as follows: time after SAA (T_SAA_HXD) > 436 s, 
the new cut-off rigidity (COR2) > 6.0 GV, elevation from 
the earth rim > 5° and that from the sun-lit earth rim 

> 20°. Non X-ray background (NXB) spectra of the 
XIS were created by sorting night earth data by COR2, 
and weighted-averaging them using the ftool "xisnxbgen" 
(Tawa et al. 2008). The estimation uncertainty in the 
NXB level is handled as a systematic error in the following 
analysis. Tawa et al. (2008) studied the error on a typical 
observation lasting for a few days. At the 90% confidence 
level, it is 6.0% and 12.5% for the sum of three XIS-FI 
(XIS-0,2,3) data and the BI (XIS-1) data, respectively. 

Another major background is the Cosmic X-ray back- 
ground (CXB). Details of the CXB estimation is shown in 
appendix 1 . Because the whole XIS region in all the three 
observations could be filled with non-thermal or thermal 
diffuse emission, we used the Lockman Hole observation 
(Suzaku observation ID, 101002010) as the CXB template. 
In addition, the northern 1/3 of the NWR pointing data 
was utilized for foreground Galactic soft thermal com- 
ponent modeling. The CXB fluctuation is estimated to 
be 11% over the XIS full FOV. Hereafter, we treat the 
background-subtracted (both the NXB and CXB) spec- 
tra as signals from Abell 3667. 

2.2.2. The XIS image 

In figure 1, we show 1-4 keV and 4-8 keV XIS-FI mo- 
saic images of the three pointings. After co-adding the 
images from XIS-0, 2 and 3, they were corrected for the 
overlapping exposure, and for vignetting effect of the mir- 
ror optics, after subtracting the NXB. Thanks to the wide 
energy range and low background of the XIS, we can see 
signals even in the NWR region which is 30' (2 Mpc) away 
from the X-ray centroid. On the other hand, no apparent 
emission associated with the radio relic is seen in neither 
band image. Detailed analysis using these images is per- 
formed in section 3.4.3. 

2.3. Data reduction of the HXD 

2.3.1. Screening and background estimation 

We processed the uncleaned event files of the HXD with 
"hxdgtigen" , and remaining events were further screened 
with HK files as follows; time after SAA (T_SAA_HXD) 

> 500 s, the COR2 > 8.0 GV and elevation from the 
earth rim > 5°. As the NXB of the PIN, we utilized 
public NXB model which is provided by the HXD team 
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Fig. 1. Mosaic XIS image of the three pointing observations, 
in (a) the soft 1-4 keV band and (b) the hard 4-8 keV band. 
See text for details. The images are smoothed by a Gaussian 
kernel with a = 0.5 arcmin. Gray-scale contours are logarith- 
mically scaled by a factor of 1.5. In the hard band image, 
the 843 MHz radio image (by SUMSS) is overlaid as a thin 
contour. Two point source candidates arc shown as open cir- 
cles in the NWR region, and the projection analysis band is 
shown as a dashed lines in panel (b); both are discussed later 
in section 3.4.3. 
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in the Suzaku website 1 . The version of the model is 
"METHOD=LCFITDT" , or "tuned". Observations of 
Abell 3667 was carried out in a period in which the lower 
threshold of the scintillator was set low, and a data trans- 
fer overflow between the HXD's analog and digital elec- 
tronics sometimes occurs. This phenomena should not 
take place in ordinary observations, and hence correspond- 
ing time regions are discarded in the default screening cri- 
teria. In fact, the effect of this phenomena in principle is 
recoverable through the dead-time correction procedure 
and thus the data should be useful. Since the HXD data 
in this paper, however, are dominated by systematic errors 
and not statistics, and most of the systematic error anal- 
ysis is performed in periods excluding these phenomena, 
we stick to the original screening criteria. Note that none 
of the results are changed significantly if we include the 
data obtained in periods with the data transfer overflow. 

In the HXD-PIN NXB documents (Fukazawa et 
al. 2008; Mizuno et al. 2008 2 ), the reproducibility of 
blank-sky observations separated into 10 ks exposures 
gave a distribution of 5.7% at the 90% confience level, 
including the statistical error of typically 3.3% or larger 
and the CXB fluctuation of 1.3%, as described later. This 
gives the NXB systematic uncertainty of 4.5% at the 90% 
confidence. In the following analysis, we utilized this value 
as the systematic error of the NXB. Because the selec- 
tion criteria of the blank-sky observations are simple and 
based only on the XIS source flux, some of them can be 
contaminated by off-center sources. In addition, 10 ks is a 
bit shorter than the exposure of our observations. Thus, 
the actual systematic error is probably smaller than this 
value. 

Detail of the CXB model used for the PIN is de- 
scribed in appendix 2. We defined the photon flux 
model as N{E) = 8.69 x 10~ 4 x (.E)" 1 - 29 x cxp(-£:/40.0) 
in photons cm~ 2 s _1 keV^ 1 FOV -1 . Here, E is the pho- 
ton energy in keV, and the model is normalized to the 
HXD-PIN opening angle (or FOV) of f^ XD = 0.32 deg 2 , 
to be combined with the HXD nominal response. The 
13.2% systematic difference of the XIS and PIN (Ishida et 
al. 2007) 3 is also accounted in the following fitting pro- 
cess. The level of the CXB fluctuation is calculated to be 
18%, which corresponds to 1.3% of the CXB+NXB in the 
10-40 keV band. 

We also checked whether the PIN data suffered from 
thermal noise in the lower energy band by comparing 
the actual earth occultation data with the correspond- 
ing model NXB spectra. Since the detector temperature 
was relatively high at —12 C° in these observations, the 
thermal noise is seen as a steep rise below 13 keV. In the 
higher energy band, the two spectra are consistent within 
statistics. 

2.3.2. Signals of the three pointings 

In figure 2, the 13-40 keV band light curves of PIN for 
a total of 4 days of observations are presented. The PIN 

1 http:/ /www. astro. isas.jaxa.jp/suzaku/analysis/hxd/pinnxb/ 

2 Suzaku-mcmo, JX-ISAS-SUZAKU-MEMO-2008-03. 

3 Suzaku-raemo, JX-ISAS-SUZAKU-MEMO-2007-11. 
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Fig. 2. The HXD light curves in the 13-40 keV range. Black 
crosses stands for the data, red for the background model and 
green for the residual signals. Time regions of the three obser- 
vations are separated by vertical dotted lines. The horizontal 
dashed line indicates the expected CXB rate. 

background varies by ~ 25%, mainly correlated with COR, 
and time elapsed after a passage of SAA (T_SAA_HXD). 
The PIN residual signal clearly exceeds the CXB level 
for all three observations. Note that the flux increases 
towards the cluster center. Small variations seen in the 
residual signal are due to uncertainties in the NXB mod- 
eling in this rather short time bin (90 minutes), and the 
NXB-subtracted count rate is consistent with being con- 
stant within each pointing. 

In Figure 3, we show raw PIN spectra from the 3 obser- 
vations, together with the NXB model spectra. Difference 
between the two spectra are also presented, and compared 
to the CXB model. In all cases, the NXB-subtracted data 
significantly exceeds the CXB. The signal is strongest in 
the center pointing, and weakest in the NWR pointing. 
The 13-40 keV PIN signal rates are (8.9 ± 1.3, 5.4 ± 1.0 
and 1.4±0.5) x 10 -2 cts s _1 in the center, the 17' offset, 
and the NWR data, respectively. The systematic error is 
1.5 x 10 -2 cts s" 1 in all pointings, when evaluated as a 
quadrature sum of those of the CXB and the NXB. Thus, 
we can claim the signal detection from Abell 3667 up to 
30 keV both in the center and 17' offset observations. 

We also estimated the point-source contribution to the 
HXD-PIN using the ROSAT-PSPC and the ASCA-GIS 
data. In short, we confirmed that the contribution from 
these sources is less than 10% and hence almost negligible. 
Details are shown in appendix 3. 

2-4- Preparation of arf files of the XIS and the PIN 

When trying to quantify the brightness and spectral 
shape of the PIN signals, we need to know the angular 
transmission (or effective area) of PIN to the correspond- 
ing hard X-ray emission. This in turn requires knowledge 
of the location and shape of the source, which is unavail- 
able in the PIN energy range. In evaluating the ICM 
contribution to the PIN band, we therefore assumed it 
to have the same spatial distribution as the ICM signals 
detected with the PSPC. Here, the ICM was assumed to 
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Fig. 3. PIN spectra of the three pointings, from (a) center, 
(b) 17' off and (c) NWR. Black crosses stand for the data, 
red for the NXB model and green for the residual signals 
(data — NXB). The expected CXB spectra are also plotted 
as green histograms. Typical systematic errors due to the 
CXB fluctuations and the NXB reproducibility are shown in 
the NWR data in blue. 

be isothermal, which is not completely true for this clus- 
ter. However, this assumption is valid as a starting point, 
since the major part of the ICM emission contributing to 
the PIN band originates from within 20' of the cluster 
center, within which region the average spectra is well ap- 
proximated by a ~ 7 keV ICM emission (see sections 3.1.1 
and 3.2). In this case, inhomogenicty in the ICM spec- 
tra is observable as a small shift in spectra of individual 
pointings. 

The PSPC image was background subtracted and cor- 
rected for the vignetting effects following the method de- 
scribed in Snowden et al. (1994). Since the ICM emission 
in the image is not observable beyond 40' from the cluster 
center, we discarded the data in this region. After nor- 
malizing the image to have a total area of 1.0 and binning 
it with a 2' x 2' mesh, we simply calculated the PIN arf file 
by convolving the point source arf of the specified mesh 
multiplied by its relative intensity. Note that the energy 
response has in principle negligible dependence on the po- 



sition within the PIN-FOV. The cross-normalization fac- 
tor of 1.132 between the PIN and the XIS is incorporated, 
as well. 

Regions beyond 26' from the cluster center, i.e. those 
as far as the XIS FOV of the NWR pointing, contribute 
only 1.3%, 4.9% and 14.4% to the normalization of arf files 
of the center, the 17' offset and the NWR pointings, re- 
spectively. Thus, 86% of the ICM photon detected in the 
NWR PIN data is expected to come from the inner region 
of the cluster. Interestingly, the normalization relative to 
the center pointing thus calculated are 0.61 and 0.14 in 
the 17' offset and the NWR pointings, respectively. This 
value matches the actually observed ratios, 0.60±0.25 and 
0.11 ±0.18, respectively. As a reference, those of a point 
source located at the cluster X-ray centroid are 0.42 and 
0.05, respectively. These results suggests that the ICM 
emission around the cluster central region must be care- 
fully handled when analyzing the NWR PIN data. 

In a similar way, we synthesized the XIS arf files with 
the ftool "xissimarfgen" , using the PSPC image as an in- 
put. Note that the input image is automatically normal- 
ized to 1.0 in the software. Thus, the spectral fitting of 
both the PIN and the XIS is normalized to the whole 
cluster, and not the FOV of each detector of individual 
observations. Since the NWR region is around the de- 
tection limit of both the PSPC and the XIS, the arf files 
for this region could be inaccurate, and they require spe- 
cial care in the following analysis. This effect is negligible 
in the center and 17' offset pointings because of the high 
X-ray surface brightness. 

To verify the accuracy of the generated arf files, we com- 
pared the XIS results to those of the PSPC. We derived 
the 0.5-2.0 keV ICM flux by integrating the PSPC sur- 
face brightness profiles given by Mohr et al. (1999) within 
40' from the X-ray centroid. The value matched within 
1% to the flux obtained by fitting the center XIS FI spec- 
tra using our arf, which is also normalized to the PSPC 
image integrated over 40'. Here, a kT = 7 keV single tem- 
perature thermal emission was assumed. Although the 
analysis in this paper is not sensitive to the Suzaku-PSPG 
cross calibration, this result supports our arf generation 
procedure. On the other hand, the XIS and PIN cross- 
calibration, to which our analysis is sensitive, is verified 
within ~ 3% accuracy using the Crab observations, in- 
cluding the mapping within the PIN-FOV. Because the 
PIN-FOV is relatively simple in its shape, no strong sys- 
tematic error is expected. Using the similar approach, the 
PIN spectra of the Abell 1060, Centaurus (Kitaguchi et 
al. 2008) and Coma clusters (Wik et al. 2008) were re- 
produced successfully by extrapolating the XIS or XMM 
spectra. From these results, we conservatively assume a 
10% uncertainty in the relative normalization of the XIS 
and the PIN arf files. 

3. Wide-band spectral analysis 

3.1. Spectral properties of the central pointing 

With the preparation made in section 2, we proceed to 
the spectral analysis of the XIS and the PIN. The central 
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Table 2. Best fit parameters of the IkT and 2kT models 
fitted to the center XIS spectra. 



IkT 



2kT 



kT\ ovf 


6.8 


a 71 +U.27 +0.07 
1 x -0.34 -0.06 






99 o+7.5 +1.1 
zz -°-5.0 -1.4 


Z* 


0.32 


QO+0.03 +0.01 
u - oo -0.03 -0.00 


N lo J 


9.5 x 10" 2 


g 4 +0.6 +0.0 1Q-2 
"•^-0.8 -0.1 * 1U 


iV hi § 




3.5±°j jjj x 10- 2 


X 2 /d.o.f. 


1139.6/1004 


1035.4/1002 



:rgy (keV) 



* Temperature of the (cooler) thermal component in the IkT (2kT) 
fit, in keV. Errors are 90% statistical and systematic, respectively, 
latter evaluated by shifting the NXB and CXB within their 
systematic errors. 

t Temperature of the hotter thermal component in the IkT fit. 
t Metal abundance. 

§ Normalization in the apec code, for the low and high temperature 
components. 



Fig. 4. The XIS-FI (black) and XIS-BI (red) center spec- 
tra, fitted with (a) lfeT and (b) 2kT thermal emission model. 
In the panel (b), two thermal components are independently 
presented as dashed histograms. 

pointing afforded the signal detection with the highest sig- 
nificance using both instruments. Their arf files are also 
the most reliable since the brightest central emission is 
well detected by the PSPC. Thus, we started the spectral 
fitting from the central pointing. We first evaluated the 
XIS and the PIN spectra individually, and then performed 
the combined fitting. 

Spectra from the three XIS-FI chips are summed up to 
form a single XIS-FI spectrum. Since the XIS-BI chip 
has a significantly different detector response, its spectra 
is separately handled. The arf files prepared in section 
2.4 are used for the XIS-FI, XIS-BI and the PIN data. 
Thus, the normalization of individual spectra should be 
the same. Practically, we fixed the relative normaliza- 
tion between the XIS-FI and the PIN data, and the 10% 
systematic uncertainty was handled separately. The nor- 
malization between the XIS-FI and BI data was set free, 
to adjust for residual mutual calibration uncertainties of 
about 5% or less. Here, the XIS-FI, the XIS-BI and the 
PIN spectra are utilized in 0.6-10.0 keV, 0.6-8.0 keV, 
and 13-40 keV, respectively. Since the data have very 
high statistics, the energy range of 1.76-1.86 keV in the 
XIS spectra were ignored to avoid uncalibrated structures 
caused by the Si edge. The XIS-FI and BI spectra are 
grouped so that each energy bin has at least 200 counts. 
The PIN spectra are grouped into fixed 12 bins, each with 
sufficient counts. 

3.1.1. The XIS 

We first applied a single temperature thermal emission 
model (hereafter IkT model) to the XIS spectra, utilizing 
the apec code in the XSPEC (vll.3.2) software package. 
Galactic absorption was modeled by the wabs code with a 
fixed column density of 4.7 x 10 20 cm -2 , as derived by the 



radio observations (Dickey & Lockman 1990). The red- 
shift was also fixed at z = 0.0556. Relative ratio of metal 
abundances are fixed to the solar value, given by Anders & 
Grevesse (1989). The temperature and the normalization 
were also set free. 

Fit with the IkT model was unsuccessful with 
X 2 /d.o.f. = 1139.6/1004, leading to null hypothesis prob- 
ability of 1.8 x 10~ 3 . The fitted temperature was about 
6.8 keV, which is generally consistent with temperatures 
given in the literature (e.g. Fukazawa ct al. 2004). Figure 
4 presents the data compared to the best fit model. In 
the IkT fitting, residuals around 1 keV and above 8 keV 
are significant. Even if we set the TVh free, the fit did not 
improve significantly, suggesting the thermal component 
modeling itself needs improvement. 

Considering the multi-temperature nature of the clus- 
ter (e.g. Briel et al. 2004), we added another thermal 
component. Metal abundances of the two apec compo- 
nents were tied to be the same, and only the temperature 
and normalization were set free (hereafter 2kT model). 
The fit improved significantly and became acceptable with 
X 2 /d.o.f. = 1035.4/1002 or null hypothesis probability of 
23%. Here, we also estimated the systematic error by 
shifting both the NXB and CXB to their highest and low- 
est possible value at the 90% confidence. The shifts ob- 
served in the fitted central value are quoted as systematic 
uncertainty. In this fitting, the hotter component temper- 
ature is suggested to be higher than 17.1 keV. Here (and 
in all of the following analysis), the errors are derived as a 
quadrature sum of the statistical and the systematic ones 
when obtaining the final result, while the two components 
are separately presented in the tables for clarity. When 
fitted using a model with a thermal and a power-law com- 
ponent (hereafter IkT + PL model), the fit also improves 
providing x 2 /d.o.f. = 1057.8/1002 or null hypothesis prob- 
ability of 10.1%. The photon index derived is hard with 
T = l.OtoJ- Before further modeling the ICM, we will 
check the HXD-PIN spectra as it should contain informa- 



No. ] 

Table 3. Spectral parameters fitted to the HXD spectra of 
the center pointing in the i3-40 keV range. 





T or kT* 


flux or norrrr 


xVd.o.f. 


PL model 
lkT model 


o c+l.U +0.5 
°- u -0.8 -0.3 
o c+5.6 +1.6 
°- J -2.9 -1.8 


q 9 +5 +6 
oz -5 -6 
7 o + 12.1 +2.2 
'•°-4.1 -1.0 


6.2/10 
6.2/10 



* Photon index V or thermal component temperature (in keV) of 
the power-law or lkT fitting, respectively. 

t The 10-40 keV flux of the power-law component, in 10 -12 erg 
s -1 cm -2 , or the normalization in apec model, for the low and high 
temperature components. 



tions about hotter or harder components. 

3.1.2. The HXD-PIN 

The signal detected by the PIN is significant compared 
to any possible systematic errors. We then studied the 
shape of the PIN spectrum in the 13-40 keV range. Here, 
we used the HXD nominal position response, and did not 
use the arf files generated at section 2.4 yet. When we 
simply fit a power-law model to the PIN spectrum, the 
resultant photon index Y is as soft as 3.6+J'g. If we try 
the lkT model with metal abundances fixed at 0.3, the 
resultant temperature is derived as 8.5I34 keV. Detailed 
results are presented in table 3. At the first glance, no 
strong signature of spectrally hard component is suggested 
in the PIN spectra. We also note that, although the error 
bar is large because of both the statistical and systematic 
errors, the photon index T is softer than those suggested 
from the radio observation, i.e. T m 2.1 and the XIS spec- 
tral fitting with lkT + PL model. On the other hand, the 
PIN flux level provides important results, as shown in the 
next section. 

3.1.3. Combined XIS- HXD analysis 

Following the spectral analysis individually performed 
to the XIS and the PIN, in this section we give the results 
of a wide-band spectroscopic analysis, using the two spec- 
tra combined (the XIS+PIN combined spectra). In order 
to connect the two instruments and derive the ICM com- 
ponent, we utilized the arf files generated in section 2.4, 
and fixed the normalization of the XIS-FI and the PIN 
spectra. 

We first applied the lkT model. Procedures and fitting 
parameters are the same to those of the XlS-only analysis. 
The fit again gave a temperature of ~ 6.9 keV. However, 
as shown in figure 5 (a), the lkT model failed to explain 
the PIN data, as well as the XIS data in the energy band 
above 7 keV. As a result, the fit is unacceptable with 
X 2 /d.o.f. = 1173.4/1016. What is more, the observed PIN 
counts are about a factor of 2 larger than the best-fit lkT 
model. This excess is larger than any possible calibration 
errors in the arf, and any systematic error associated with 
the NXB and CXB, as well. 

To characterize the hard excess, we next applied the 
2kT model. Relative normalization of the two compo- 
nents are the same in the XIS and the HXD-PIN spec- 
tra. In other words, we assume a simple model in which 
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Fig. 5. Wide-band spectra of the center pointing, fitted with 
the (a) lkT, (b) 2kT and (c) lkT + PL model. In the panel 
(a), typical 90% confidence systematic error in the NXB and 
CXB are also shown for reference, in cyan, blue and magenta 
corresponding to the XIS-FI, BI and PIN, respectively. Panel 
(a) is shown together with the residual ratio, while the other 
two are with the delta-chi distribution. 

the ICM emission uniformly contains the two compo- 
nents. Since 75% of the ICM emission accumulated in 
the center PIN data comes from within 12.5' from the 
X-ray centroid, which is mostly covered by the center 
XIS FOV, this model is valid as a simple approximation. 
As shown in figure 5 (b) and table 4, a good fit with 
X 2 /d.o.f. = 1046.5/1014 or a null- hypothesis probability of 
23.3% was obtained. Interestingly, both the XIS and the 
PIN hard excess consistently disappeared. Considering 
both the statistical and systematic errors (table 4), this 
hard component, if interpreted as thermal, is character- 
ized by a temperature of kT^i ranging from 13.2-25.8 keV. 

As another candidate, we fitted the lkT + PL model 
to the XIS+PIN combined spectra. The fit was improved 
compared to the lkT model, and marginally acceptable 
giving 4.8% in null hypothesis probability. As shown in 
table 4, the photon index was hard with r = 1.4 and the 
inferred power-law flux in 10-40 keV was ~4x 10~ n erg 
s" 1 cm -2 , normalized to the cluster as a whole. 

The 2kT model is the simplest one to evaluate the 
multi-temperature nature of the cluster ICM. Then we 
also tried a model with three apec components (hereafter 
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Table 4. The best fit parameters of the 2kT and lkT + PL 
models fitted to the XIS+PIN combined spectra of the center 
pointing. Columns are similar to those shown in table 2. The 
systematic error includes the XIS-PIN arf uncertainty of 10%, 
in addition to those of the CXB and NXB. 



2kT lkT+PL 



kT\ ow 


d eq+u.yy +o.oi 

*-O»_Q.40 -0.08 


U.Ud_ 20 -0.28 


kT hi 


iq 9+ 4 - 7 + 4 - 6 

13 - z -4.0 -4.5 




Z 


n on+0.03 +0.00 
u - oy -0.03 -0.01 


n or+0.03 +0.01 
u - oo -0.03 -0.01 


N lo 


o -,+0.6 +0.2 ln -2 
D - L -0.8 -0.4 x 1U 


o O+0.5 +0.5 v in-2 
°- z _0.4 -0.3 x 1U 


N hi 


o fi +0.8 +0.4 1n -2 
t3 -°-0.6 -0.0 x 1U 




TpL* 




1 QQ+0.10 +0.00 
l.O»_ !7 _o.02 


iW 




o 7 +7 +14 
°'-6 -13 


X 2 /d.o.f. 


1046.5/1014 


1089.8/1014 



* Photon index of the power-law component. 

' 10-40 keV flux of the power-law component, in 10 — 12 erg s — 1 

_9 

cm . 



Table 5. The best fit parameters of the 3kT model fitted to 
the XIS+PIN combined spectra of the center pointing. 



3fcT 



fellow 


o co+0.56 +0.00 
°- uo -1.14 -0.86 






7.5 (fixed) 




fcT hi 


97 o+oo +oc 
z ' • J -8.2 -0.0 




Z 


n qo+o.02 +0.00 

u -°°-0.02 -0.02 




N lo 


3.0±i-l ±Hxl0" 


2 




4-7±i;i x 10- 


2 




2.i±i-i x io- 


2 


X 2 /d.o.f. 


1046.5/1014 





* A fixed temperature of the medium hot component. 

3fcT model). In this case, we fixed one of the three tem- 
peratures to 7.5 keV to make the fit converge. The value 
refers to the lkT fit to the 17' offset data shown in the 
next section, and is also consistent with one of the major 
components seen in the XMM mosaic (Briel et al. 2004). 
As shown in table 5, the fit was acceptable with a null 
hypothesis probability of 27%. A hot component with 
a temperature fcThi = 27.3^^° 2 keV is required, but with 
~ 40% smaller normalization compared to those of the 
hotter component of the 2fcT model. If we shift the tem- 
perature ranging from 6.5 to 8.0 keV, still the fit is ac- 
ceptable and the parameter value is generally similar. In 
other words, we still need a component as hot as ~ 20 
keV, which could also be replaced by a V ~ 1.4 power law 
component. 

The temperature map drived from XMM mosaic shows 
several hot regions with fcT > 8 keV around a few minutes 
north-west to the X-ray centroid (see figure 2 of Briel ct 
al. 2004). Using Chandra, the signature of a component 
as hot as > 10 keV is also suggested in the same region 
(Vikhlinin et al., private communication). Due to pro- 
jection of emission components, these value should be re- 
garded as a lower limit, and thus a very hot (> 13.2 keV) 
component apparently is required in the cluster. In the 
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Fig. 6. Wide-band spectra of (a) the 17' offset pointing fitted 
with lfcT model, (b) plotted with 2kT best-fit model of the 
central pointing, (c) The NWR spectra with the lfcT best- 
fit model to the XIS. Cyan, blue and magenta plots show 
the typical systematic error of the XIS-FIs, BI and the PIN 
spectra. See text for details. 

following analyses, we regard the 2fcT model as an upper 
limit of the ICM contribution to the PIN-band, which is 
the most important issue in this paper. Using additional 
long-exposure XMM data to be obtained soon, we will 
further pursue this issue in the next paper. 

Through these analyses, we conclude that the excess 
hard component suggested in the XIS data is clearly re- 
quired in the PIN spectra as well. The emission is likely to 
be of thermal origin, a very hot component with tempera- 
ture exceeding 13.2 keV, or less likely a hard (r ~ 1.5) 
power-law component. We also checked the Fe-K line 
profile, and found that it was consistent with either the 
lkT + PL, 2kT and 3fcT interpretations. In any case, the 
analysis of the other observations, i.e. the 17' offset and 
the NWR pointings, will be helpful to identify the origin 
of the hard emission. The XIS and PIN data in these 
pointings contains spatial information of the excess com- 
ponent. What is more, the strong north-west radio relic is 
observed with the deep (81 ks in the XIS) NWR pointing. 

3.2. Spectral properties of the offset pointings 

When analyzing the offset pointings, we have to take 
into account that, although the XIS spectra are of local 



No. ] 

Table 6. Similar to table 4, but for the 17' off and the NWR 
offset spectra fitted with the IkT model. 



17' off* nwr; 



kT [keV] 


7 7-1+0.33 +0.38 
' • ' x -0.32 -0.47 


o R+0.4 +0.5 
°-"-0.4 -0.6 


Z 


n n fi +0.06 +0.01 
u - zu -0.06 -0.01 


n ifi+0.16 +0.01 
u - lu -0.14 -0.00 




9.3t o; 2 2^: 2 xio- 2 


7.4t° o: f x io- 2 


r7d.o.f. 


234.1/220 


120.5/89 



* Fitted to the XIS+PIN combined spectra, 
t Fitted only to the XIS spectra. 

t Normalization in apec model. Note that it is normalized to the 
whole cluster, not within the observation FOV, using the arf files 
generated from the PSPC image. 



origin (from within the FOV), the PIN spectra is largely 
contaminated from the emission from the cluster center. 
For example, 59% and 37% of the ICM emission accu- 
mulated in the 17'offset and the NWR pointings, respec- 
tively, comes from within 12.5' from the X-ray centroid, 
i.e. around the center. Therefore the XIS and the PIN 
spectra wil not necessarily show the same results, although 
they should be roughly consistent with one another as a 
whole. 

The XIS spectra of the 17' offset pointing can be rep- 
resented using the IkT model, with kT = 7.6^4 keV and 
x2/d.o.f. = 221.9/208, or a null hypothesis probability of 
24.3%. Note that we utilized the arf files generated in sec- 
tion 2.4 so that the (energy dependent) vignetting effect 
is fully taken into account. For the same reason, the fitted 
flux is normalized to the reference image used in the arf 
generation (in this case the whole cluster), not just the 
flux within the FOV. When we use the XIS+PIN com- 
bined spectra, the IkT model is acceptable as well, with 
a null-hypothesis probability of 24.5%. Results of the fit- 
ting are presented in table 6. As shown in figure 6 (a), 
however, the model slightly (by 20%) under-predicts the 
data in the PIN band. The residual resides right on the 
90% upper-limit of the CXB and NXB systematic error. 

The 17' offset data is also consistent with the 2kT pic- 
ture of the center pointing. In fact, when we applied the 
2kT best fit model of the center pointing to the XIS-PIN 
spectra, a marginally acceptable fit was obtained by shift- 
ing its global normalization only by —3% (see also panel 
(b) of figure 6). The fit gives x2/d.o.f. = 260.8/222 or 
a null-hypothesis probability of 3.8%. In this case, the 
model slightly (by 20%) over-predicts the PIN data. One 
possibility is that the cluster as a whole can be represented 
by the 2kT model, if one includes the systematic errors. 
The result can also be understood if the very hot compo- 
nent resides within the center XIS FOV, and the majority 
of the 17' offset XIS FOV is covered by a single - 7.5 keV 
thermal component. 

In both the XIS and PIN spectra of the NWR point- 
ing, the CXB and NXB systematic errors are not negligi- 
ble, as shown in the panel (c) of figure 6. Here, we used 
the southern 1 /3 region of the NWR pointing in the XIS 
analysis (hereafter 1/3 NW-south region) because the sys- 
tematic errors make spectral analysis difficult in the other 
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Table 7. Spectral parameters of the power-law model fitted 
to the HXD spectra of the 17' offset and the NWR pointings 
in the 13—40 keV range. 



Pointings 


r 


flux* 


x'Vd.o.f. 


17' offset 
NWR 


o A +L:2 +1.0 
°-^-1.0 -0.5 
o 4 +2.0 N/A t 
°'^-1.5 -0.9 


on+4 +6 

r 1+1.7 +6.3 
°-- L -1.7 -5.1 


8.4/10 
7.5/10 



* 10-40 keV flux of the power-law component, in 10 -12 erg s — 1 
cm -2 . 

t N/A means the value was not available. In this case, the lack of 
photon made it impossible to restrict the spectral shape. 



regions. We restrict our fitting to the energy band from 
0.7-7.5 keV and 0.7-6.0 keV in the XIS-FI and BI spectra, 
respectively. We first applied the IkT model to the XIS 
spectra. The best-fit temperature thus obtained is as low 
as 3.6lg;7 keV (see table 6), which is significantly cooler 
than those of the center and 17' offset data. The fit is 
marginally acceptable with a null-hypothesis probability 
of 1.5%. Noticeably large residuals appears around 0.9 
keV suggesting an additional ~ 1 keV thermal emission 
component. The XIS data around the north-west radio 
relic is further analyzed in section 3.4.3. 

The NWR PIN signal could be explained by the CXB 
and NXB systematic errors, since the level is right on their 
90% confidence upper- limit (see the panel (c) of figure 6) . 
Another reasonable possibility is the ICM contribution 
from the cluster central region to the PIN spectra. Based 
on the arf files derived in section 2.4, 37% and 49% of 
the ICM contribution to the NWR-PIN data is originated 
from within 12.5' and 26' from the cluster X-ray centroid, 
respectively, which roughly correspond to the area of the 
center and 17' offset XIS FOVs. As shown later in section 
3.4.1, the ICM contribution estimated using the XIS can 
explain all of the PIN signals in the NWR pointing. As a 
cross check, we also introduce the XMM mosaic data (e.g. 
Bricl ct al. 2004) covering the entire PIN FOV later in 
section 3.4.2, and come to the same conclusion. 

In table 7, we present the PIN spectral properties sep- 
arately fitted to the 17' offset and the NWR data. Using 
the power-law model, the photon index T becomes as soft 
as 3.5. Although the error is large, the best-fit value is 
consistent with those obtained from the center PIN spec- 
tra (sec table 3), while the flux itself is much brighter 
there. These results strongly suggest that the same emis- 
sion dominant in the center-PIN spectra is contaminating 
both the 17' offset and the NWR PIN data. 

3.3. Joint fit to the wide-band spectra of the three point- 
ings 

The PIN data of the three pointings are all mixed 
to some extent, mostly contaminated from the emission 
around the cluster center. In this section, we combine 
them to obtain the cluster overall properties. 

Analysis of the three pointings suggests that all the data 
are consistent with the simple two component picture sug- 
gested in the center pointing data, except for the XIS data 
in the NWR pointing showing rather low ICM tempera- 
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Table 8. Best fit parameters of the 2kT model and IkT+PL 
model fitted to the all-combined spectra. Columns are similar 
to those presented in table 4. 



2kT IkT+PL 



kT lo 


A 79 +().3H +0.01 
*-' z -0.60 -2.81 


°- w -0.16 -0.43 


kT hi 


-, A 7 +3.1 +10.4 
—3 2 — 7 2 




Z 


n Qc+0.02 +0.00 

U.OO_q q 2 _o.04 


n cm +0.02 +0.02 
u -°^-0.02 -0.02 


N lo 


5.6+°i +0 5 8 x 10- 2 


o r+0.4 +0.8 1f) -2 
°-°-0.3 -0.8 x 1U 


N hi 


+ 1.3 +4.8 1n _2 




Tpl 


1 AR+ - 10 + 001 
1 - w -0.16 -0.02 


-Fpl 




1+5.3 +22.5 
z0 " L -4.8 -21.7 


X 2 /d.o.f. 


1320.4/1248 


1364.0/1248 



ture. Because of its low surface brightness, thermal con- 
tribution from this region is minor in the PIN band. Thus, 
we conducted a spectral analysis to the data combining in 
total 5 data sets (hereafter all-combined data), i.e. three 
PIN spectra and the XIS spectra of the center and 17' 
offset pointings. Following the strategy employed at sec- 
tion 3.1, we fixed the relative normalization between the 
center XIS-FI spectra and all the PIN spectra, while those 
between the XISs (offset FI data and all the BI data) were 
set free to handle small (< 5%) residual uncertainties. As 
noted, the 10% systematic error between the XIS-FI and 
the PIN data were handled separately. 

As was the case with the center pointing, the lkT model 
fitting to the all-combined spectra was unsuccessful, giv- 
ing x 2 /d.o.f. = 1430.7/1250 or a null-hypothesis probabil- 
ity of 3 x 10 -4 . Thus, we applied the 2kT model to the 
all-combined spectra. 

The 2kT fit was almost acceptable as shown in figure 
7 and table 8. The null- hypothesis probability of the fit 
was 7.6%. Partly because of the slight discrepancy of the 
spectral shape around 1 keV between the center and the 
17' offset poitings, the hotter component temperature was 
not determined well when the NXB, the CXB and the arf 
systematic errors were included in the fit. However, the 
best-fit value of the hotter component temperature kT\^ is 
still 14.7 keV, which is consistent with the results from the 
center pointing data (see section 3.1.3). The 3fcT model 
gives acceptable fit, as well. 

The IkT + PL model fitting to the all-combined spec- 
tra gave results similar to those of the center pointing 
(table 8). The fit itself was marginally acceptable with 
a null-hypothesis probability of 1.2%, significantly worse 
than that of the 2kT model fitting. Although the poor fit 
could be due to an overly simplified model, the fact that 
the HXD spectra have a shape as soft as T = 3.5 strongly 
indicates emission of thermal origin. Thus here we con- 
clude that the all-combined fitting also prefers the 2kT 
(or 3fcT) picture. Another picture that the very hot com- 
ponent resides only around the cluster center may apply 
as well, considering the systematic errors. In any case, 
the NWR-PIN data is well explained by the contribution 
from the cluster central portions, suggesting there is no 
strong signal of the IC emission in the data. 




energy (keV) 

Fig. 7. The all-combined spectra fitted with the 2kT model. 
For simplicity, the XIS-BI spectra are not shown, although 
used in the actual fitting. The XIS and the PIN spectra of 
the center pointing are shown with black crosses, and those of 
the 17' offset pointing in red. Green crosses give the PIN spec- 
trum of the NWR pointing. Two components of the model 
is presented individually as a dashed line only for the center 
data. 

3-4- Search for the IC signal from the north-west radio 
relic 

The north-west radio relic region is mostly covered by 
the deep (81 ks in the XIS and 48 ks in the HXD) NWR 
pointing. Although there is no strong evidence for the IC 
emission in the analysis so far, there should be some X-ray 
and hard X-ray emission from the relic. In this section, 
we focused mainly on this region to detect or determine 
the upper-limit on the IC emission. 

3.4-1- The HXD upper limit combined with the XIS re- 
sults 

We evaluated the flux of the possible non-thermal hard 
X-ray emission using the 13-40 keV PIN signal of the 
NWR pointing. In table 9, signal rates are summarized 
with their errors. To estimate the ICM components, we 
used both the 2kT model fitted to the center XIS-PIN 
combined spectra and the lkT model fitted to the 17' 
offset XlS-only spectra. The 17' offset PIN data was not 
used, because it may contain some signal from the north- 
west radio relic. 

Using the center-2fcT model, we converted the residual 
signal in the NWR-PIN data by assuming a T = 2.1 fixed 
power-law emission located at the FOV center. The result 
is (-5.1±1.8±5.8) x 10" 12 erg s _1 cm" 2 , i.e. < 1.0 x 10~ 12 
erg s^ 1 cm -2 in the 10-40 keV band. Here again, the CXB 
and the NXB systematic errors and the 10% systematic 
error in the arf files are taken into account. Note that the 
arf files are used to estimate the ICM contribution, but not 
in the hard excess flux evaluation. Using the lkT model 
fitted to the 17' offset XIS spectra, the power-law flux is 
estimated to be (-0.7 ± 1.8 ±5.7) x 10~ 12 erg s" 1 cm" 2 , 
i.e. < 5.3 x 10~ 12 erg s -1 cm~ 2 in the 10-40 keV band. 
Since the contribution from the central portion and 17' 
offset region are both significant, actual upper-limit on the 
power-law emission will be somewhere between these two 
values. Here, we recognize the latter results, < 5.3 x 10~ 12 
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Table 9. Event counts in the 13-40 keV NWR-PIN data, 
compared to various error components. 



Components 


counts s 1 


Signal 


(14.1 ±4.9) x 10~ 3 


4.5% of NXB* 


±15.0 x 10~ 3 


18% of CXBt 


±4.8 x 1CT 3 


ICM 2feT * 


(28.4 ±3.0 ±2.8) x 10~ 3 


ICM lfeT t 


(16.0 ±2.0 ±1.6) x 10~ 3 


lCM X MM t 


(14.7 ±1.5 ±1.5) x 10" 3 


Signal 2 fcT S 


(-14.3 ±4.9 ±16.3) x 10~ 3 


SignalifeT § 


(-1.9±4.9±16.0) x 10~ 3 


SignaljfMM § 


(-0.6±4.9±15.9) x 10" 3 


FlkT PL # 


-5.1±1.8±5.8 


F\kT PL* 


-0.7±1.8±5.7 


Fx MM PL* 


-0.2±1.8±5.7 



* Estimated NXB systematic error at the 90% confidence. 
' Estimated CXB fluctuations at the 90% confidence. 

i ICM emission component estimated using the 2kT model fitted to 
the center XIS-PIN combined data, the IkT model fitted to the 17' 
offset XIS data, and to the XMM weighted data. Values are shown 
with 90% fitting error including the NXB+CXB systematic effect in 
quadrature sum and 10% arf calibration systematic error. 
§ Non-thermal signals in counts s _1 with 90% statistical and sys- 
tematic errors. The latter value is the quadrature sum of those 
associated with NXB, CXB and the ICM. 

# The 10-40 keV band flux in 10 -12 erg s _1 cm -2 derived from 
the "Signal" column, converted to a T = 2.1 power law emission. 

Table 10. Best fit parameters of the IkT model fitted to 
the XMM-pn spectra, weighted with the HXD-PIN vignetting 
function of the NWR pointing. 



kT* 


Abun.t 


xVd.o.f. 


' - Jt3 -0.42 


n 99+0.04 


116.6/123 



* Temperature of the apec model. 

* Metal abundance of the apec model. 

erg s _1 cm -2 , as the conservative upper-limit flux of the 
excess hard X-rays from the north-west relic. 

3.4-2. The HXD upper limit combined with the XMM 
results 

Mapping observations performed by XMM, as described 
in Briel et al. (2004), provides independent information of 
the ICM and its temperature distribution. With its larger 
FOV, the XMM data has a merit on observing almost 
all of the ICM emission covered by the HXD-PIN FOV, 
while the XIS has the merit of low background and simpler 
cross-calibration with the HXD. 

The XMM spectra, background and arf files for this 
analysis are generated using the method described in Wik 
et al. (2008). In summary, we defined 10 regions with 
nearly constant PIN effective area, separated by curves 
with 0%, 10%, 20%, ... 90% of the value for the central 
point source, as shown in figure 8. Then we extracted 
XMM EPIC-pn spectra from these region in the standard 




16:00 14:00 12:00 20:10:00 08:00 06:00 
Right ascension 

Fig. 8. The wavelet smoothed XMM mosaic image of Abcll 
3667 with a logarthmic color scale (Briel et al. 2004). The 
white square is the total FOV of the PIN for the NWR ob- 
servation. The interior white regions are curves of constant 
PIN effective area, ranging from 10% to 90% of the effective 
area for a central point source. The cyan contours are from 
the SUMSS 842 MHz radio image of the NW radio relic. Sec 
text for detail. 
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Fig. 9. The NWR-PIN spectra combined with the XMM 
spectra. The best-fit lfcT model is shown. 

manner, using the XMM-mosaic data and adopting the 
XMMSAS v7.1 calibration. The background spectrum is 
derived from the datasets compiled by Read & Ponman 
(2003), and we also apply consistent flare cleaning criteria 
to both the data and the background. The arf files used 
for the XMM data include the corrections for exposure 
and vignetting for the XMM EPIC-pn. 

By adding the 10 signal spectra multiplied by the PIN 
vignetting function, and correcting the XMM-XIS cross 
normalization of 15%, we obtain the XMM spectrum that 
would be observed with XMM but if it had the spatial 
spatial response as the Suzaku PIN detector. The nor- 
malization of the spectrum is based on the response of 
the PIN to a point source at the center of the PIN field 
of view; thus, the fluxes are correct if the PIN spectra 
are fit using the arf appropriate for a central point source. 
To avoid the contribution of the soft emission components 
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Fig. 10. The XIS 1-4 keV mosaic image (gray scale) with 
843 MHz image overlaid (contours). Regions 2' around the 
two point sources are shown as solid circles. The north-west 
relic region, used in the spectral analysis, is shown with thick 
dashed line. 
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Fig. 11. The XIS spectra of the north-west relic region fit- 
ted with 2kT thermal model. Note the cool component at 
~ 0.9 keV. 

to the results, we performed the spectral fitting above 2 
keV, and also masked the energy-region in which the back- 
ground lines (mostly from Ni) are prominent. 

In figure 9, we present the XMM-pn spectra fitted with 
the lkT model. Parameters are summarized in table 10. 
The NWR-PIN spectra is remarkably consistent with the 
XMM-best-fit model superposed, as also shown in the fig- 
ure. Thus, using the XMM, we independently obtained 
a result that no additional hard X-ray component is re- 
quired in the NWR-PIN data. As shown in table 9, the 
power-law flux is estimated to be (— 0.2±1.8±5.7) x 10" 12 
erg s _1 cm -2 , i.e. < 5.8 x 10~ 12 erg s" 1 cm" 2 in the 10-40 
keV band. 



3.4-3. The XIS image and spectral analysis 

Most of the north-west radio relic is covered by the XIS 
FOVs of the 17' offset and the NWR pointings. Since 
the background of the XIS is low and stable, the data 
contains good informations on the IC emission from the 
relic itself. We then analyzed the XIS spectra obtained 
from the region as shown in figure 10, which covers 90% 
of the 847 MHz radio flux (hereafter the north-west relic 
region). We also examined the projected brightness profile 
in three energy bands, i.e. 1-2 keV, 2-4 keV and 4-8 keV, 
to search for any signal of IC emission. 

In the north-west relic region and areas further north, 
the diffuse X-ray emission is weak and dimmer point 
sources can be identified easily. As shown in figure 10, 
we identified two sources with flux 1.5 ~ 3.5 x 10" 14 erg 
s -1 cm -2 in the 2-10 keV band, and masked out 2' 
around them. When wc perform the same operation to 
the Lockman Hole observation, we found two sources as 
well. There the overall CXB level decreases by ~ 10% af- 
ter masking these sources. This shift is fed back into the 
CXB estimation of the NWR analysis here. Since the up- 
per cutoff flux is changed to S c ~ 2 x 10" 14 erg s" 1 cm" 2 
from ~lx 10" 13 erg s" 1 cm" 2 , the CXB fluctuation level 
is reduced to 15.9% from the original 21% (see also ap- 
pendix 1). 

In figure 11, the XIS-FI and BI spectra of the north-west 
relic region data obtained after masking the sources are 
shown. Both the NWR and 17' offset data were added. To 
estimate the flux associated with the relic, we calculated 
the detector response to a widely distributed flat emission 
accumulated within the north-west relic region, and then 
scaled it to the sky coverage (f2 e = 0.04 deg 2 ) of the re- 
gion. When fitted with a lkT model, the fit is marginally 
acceptable, giving x 2 /d.o.f. = 158.1/127 or 3.2% in null 
hypothesis probability. The temperature is 5.8 ±0.7 keV, 
between the two values obtained in the XIS data of the 
two pointings (see section 3.2). To handle the significant 
hump at around 0.9 keV, we added another thermal com- 
ponent (the 2kT model). Then we obtained almost ac- 
ceptable fit with x 2 /d.o.f. = 145.9/125 or 9.7% in null 
hypothesis probability. The hot component temperature 
is derived as kT^ — 7.0t 2 'g keV and that of the cool one 
as kT\ = 1-03^0 92 keV, while the common abundance is 
0.17lg 12 solar. The significance of the Fe-K lines is weak, 
which is not surprising considering the low surface bright- 
ness. Since the flux of the cool component at ~ 0.9 keV is 
similar to that of the CXB model, the emission could be 
either a slight hump in the "foreground Galactic compo- 
nent" or some cool component really associated with the 
cluster periphery 

Since the region coincides in position with the diffuse ra- 
dio emission, the detected X-rays might be of non-thermal 
origin. As an extreme case, we fit the spectra with a 
power-law model. Because of the 0.9 keV hump, which 
indicates a contribution from kT ~ 1 keV thermal compo- 
nent, we utilized the energy band above 1.5 keV in this 
fitting. We also limit it to below 6.0 keV since the NXB 
systematic error dominates above this energy. The results 
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Table 11. Spectral parameters of the power-law model fitted 
to the XIS spectra of the north-west relic region . 



flux* 



FD7T 
-0.1 



-o.i 



0.62 



A). VI 



xVd.o.f. 



0.10 -0.19 



^ 65.1/66 



10-40 keV flux of the power-law component, in 10 



are listed in table 11. The power-law fit is acceptable 
with a null- hypothesis probability of 51%. Taking into ac- 
count the large systematic errors associated with the NXB 
and CXB (see figure 6c), wc cannot distinguish whether 
the emission is thermal or non-thermal from spectroscopy 
only. 

In figure 12, wc present the projected X-ray count-rate 
map in the energy bands of 1-2 keV, 2-4 keV and 4-8 
keV derived from all the three XIS-FI detectors. Here, the 
BI chip is not utilized because of its higher background. 
Regions around the two point source in the NWR data 
are excluded in the plot. Overlaid is the profile of the 
SUMSS 843 MHz radio image, which clearly shows the 
relic, located ~ 30' (2 Mpc) away from the X-ray centroid. 
From the 1-2 keV and 2-4 keV profiles, it is clear that 
the cluster emission was detected out to the region 40' 
from the cluster X-ray centroid, which corresponds to 2.6 
Mpc in projected distance. On the other hand, there is 
apparently no excess signal associated with the radio relic, 
smoothly connecting the inner ICM X-rays to the CXB 
component at the outermost region. Thus, a significant 
fraction of the X-rays in the NWR region was naturally 
assumed to be of thermal origin. 

Using this plot, we roughly evaluated the upper- limit 
of the projected non-thermal count rates associated with 
the relic. Here, wc assumed that the IC emission has the 
same morphology to that of the radio emission, and the 
2-8 keV X-ray counts within region B in figure 12 are the 
upper-limit to the IC emission. The CXB is estimated 
from region C, and its fluctuation is calculated to be 22% 
on both the region A and B. Taking the statistical error 
and the CXB fluctuation into account, relative count-rate 
of region B to that of region A is derived as 0.48^q 17 
while that of the radio flux is 1.17. Thus, we concluded 
that 40^16^1 i- e - l ar g er than 24%, of the emission in 
region A+B is of thermal origin. Since the region A+B 
coincides with the norht-west relic region, we combine this 
result with the power-law fitting results (table 11). In 
addition, because the radio emission within the spectral 
fitting region only covers 90% of the total relic flux in the 
SUMSS image, we multiply the XIS flux by a factor of 1.1 
to account for any missing nonthermal flux. Finally, we 
obtain the 10-40 keV upper limit IC flux for the north- 
west radio relic whole as 7.3 x 10 13 erg s 1 cm 
allowing T to vary, at the 90% confidence level. With this 
flux, the photon index is required to be T = 1.6. If we 
fix it to r = 2.1, the 10-40 keV upper limit flux become 




3.9 x 10" 13 



erg s 1 cm 2 . 



20 30 
arcmin 

Fig. 12. Projected X-ray image in the energy band of 1-2 
keV, 2—4 keV and 4-8 keV in red, green and blue crosses, 
respectively. The width of the projection region is 12' along 
the merger axis, as shown in figure 1. Black histogram stands 
for SUMSS 843 MHz radio image obtained via Skyview service 
by NASA, projected in the same manner but with arbitral 
normalization. The region A, B and C are those used for IC 
upper-limit analysis. See text for detail. 

4. Discussion 

4--1- Short summary of the analysis results 

We analyzed the Suzaku mapping observation data of 
a merging cluster Abell 3667. The XIS detects emission 
from the cluster out to ~ 40' or 2.6 Mpc from the cluster 
X-ray centroid. The PIN spectra of the three pointings 
commonly prefer soft-shaped spectra, suggesting that the 
thermal contribution from the cluster ICM is the major 
origin of these signals. From the XIS+PIN combined anal- 
ysis of the center pointing, a multi-phase picture such as 
the 2kT model, the lkT + PL model or the 3kT model is 
suggested. In the former case, fclhi is derived to exceed 
13.2 keV, with the best fit value around 20 keV, while kT\ 
is around 4.7 keV. Using the 3kT model with one of the 
temperature fixed at 7.5 keV, the fit is also acceptable, 
while a component wiht fcThi ~ 20 keV is still required. 
This result suggests the existence of a very hot (kT ~ 20 
keV) component around the cluster center. 

The XIS spectra of the 17' offset pointing can be fit- 
ted using the lkT model with kT ~ 7.5 keV, while that 
of the southern one third of the NWR pointing (the 1/3 
NW-south region) requires a cooler value, around 4 keV. 
The PIN signal in the NWR pointing, as well as the 17' 
offset pointing, can be explained by the ICM contribu- 
tion with either the 2kT picture suggested by the center 
pointing or the lkT picture fitted to the 17' offset XIS 
data. The shape of the PIN spectra of the center and 17' 
offset pointings are too soft to account for the T = 2.1 IC 
emission suggested from the radio observation. Thus, the 
wide-band Suzaku X-ray spectra do not show any signa- 
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ture of the IC emission associated with the strong (3.7 Jy 
at 1.4 GHz) radio synchrotron emission. 

Using the HXD, the overall upper-limit 10-40 keV 
flux on the non-thermal emission within a 34' x 34' FOV 
around the north-west radio relic is derived to be 5.3 x 
10" 12 erg s _1 cm " 2 with the lkT ICM modeling. The 
ICM modeling using the XMM-pn spectra gives a simi- 
lar upper limit. Using the XIS, the upper-limit on the 
IC emission from the radio relic itself is derived to be 
3.9 x 10 -13 erg s _1 cm ~ 2 , when converted into the 10- 
40 keV band and assuming r — 2.1. If we allow the in- 
dex to be free, i.e. assume that there is a spectral break 
somewhere above 10 keV, the upper-limit value became 
7.3 x 10 -13 erg s _1 cm -2 . 

4-2. The magnetic field in the north-west radio relic 

By combining the radio and X-ray observations, here we 
estimate the magnetic field strength in the north-west ra- 
dio relic. From the radio observation, the relic has a flux of 
3.7 Jy at 1.4 GHz (Johnston-Hollit 2008). Its radio fluxes 
at 85 MHz, 843 MHz, 1.4 GHz and 2.3 GHz are reported 
to roughly follow the power law with T = 2.1 (Roettgering 
et al. 1997). With a magnetic field of 1 fj,G, the 85 MHz 
radio emission requires electrons with a Lorentz factor 7 
of ~ 5 x 10 3 , which in turn emits IC photon at around 
20 keV. In the case of 1.4 GHz radio emission, they are 
7 ~ 2 x 10 4 and 350 keV, respectively Thus, we will as- 
sume that the photon index is similar to r = 2.1 above 20 
keV, when the magnetic field is stronger than 1 /xG. 

The radio flux at the 80-320 MHz band, which cor- 
responds to the ~ 20 keV IC emission with B = 1 /xG, 

is derived as F 8 q™ 2 o MHz = 1-8 x 10 -13 erg s" 1 cm -2 , 
assuming a T = 2.1 power- law emission and a flux of 
3.7 Jy at 1.4 GHz. The XIS upper-limit right on the 
north-west radio relic region is derived as -F10-40 kcV < 
3.9 x 10 -13 erg s _1 cm~ 2 , again assuming a simple T = 2.1 
power law. By dividing these two values, we obtain 

F Sync ^ tj b 
F lC ~ C/CMB 

_ B 2 1.8 xl0~ 13 
-(^cmb) > 3.9 x 10- 13 ' 

Here, -B^ MB = 3.2 ixG is the equivalent magnetic field of 
the energy density of the CMB. Thus, the magnetic field 
in the north-west radio relic is inferred to be B > 2.2 /xG. 

Since the radio emission is detected only down to 85 
MHz, the IC spectra below ~ 20 keV might have a break. 
For the electrons with 7 = 5 x 10 3 within a thermal gas den- 
sity of ~ 1.3 x 10~ 4 cm~ 3 (derived from the XIS results), 
the dominant cooling processes are the IC and the syn- 
chrotron emissions. These give a cooling time of ~ 4 x 10 8 
years. If the duration of electron injection is shorter than 
this time-scale, the electron distribution below this en- 
ergy is not affected by the cooling, and the resultant IC 
(and synchrotron) photon index can be harder, typically 
1.5 < r < 2.1. From the T free power-law fitting to the 
XIS spectra obtained from the north-west relic region, wc 
know that the upper-limit IC flux in the 10-40 keV band 



is 7.3 x 10~ 13 erg s _1 cm~ 2 . In this model, the photon 
index is 1.6 around 4 keV, and breaks to 2.1 above ~ 20 
keV. In this case, the magnetic field strength is inferred to 
be B > 1.6 fiG. These results are consistent with the radio 
rotation measure observations of this relic, suggesting 3-5 
/xG (Johnston-Hollitt 2004). 

Based on the observable parameters, we calculate the 
energy density of individual components in the north-west 
radio relic. The thermal particle density and temperature 
give an energy density Uth of 1.0 eV cm~ 3 . Here, we 
simply assume that all X-rays in the north-west relic are 
of thermal origin. With a magnetic field strength of > 2.2 
/xG and > 1.6 /xG, the energy density Ub is derived as 
> 0.12 eV cm~ 3 or > 0.07 eV cm~ 3 , respectively. Here, 
on the contrary, we assume that most of the X-rays are 
of IC origin to give an upper limit on U e . If a power-law 
distribution of relativistic electrons are assumed within 

5 x 10 2 < 7 < 4 x 10 4 , the energy density of electrons U e is 
< 0.15 eV cm -3 . Here the minimum 7 is assumed to be 
the critical Lorentz factor below where ionization losses 
dominate (Sarazin 1999). Although U e strongly depends 
on the unknown electron distribution below 7~5x 10 3 , 
the calculated upper-limit value is not so far from the 
cqui-partition condition. 

The ratio between the non-thermal (Ub and U e ) and 
thermal (t/th) energy density is inferred to be as large as 
20%. Although U e has a large uncertainty, the lower limit 
on Ub is much more robust. In addition, Us/Uth > 0.07 
is a strong limit. If wc assume, for example, that about 
a half of the X-rays in the north-west relic are of thermal 
origin and the latter half is of IC origin, Uth is decreased 
by a factor of \/2 and the lower limit on the Ub is doubled. 
In this case, the ratio becomes as high as 19%, excluding 
the U e contribution. Our results arc generally consistent 
with those by Rottgering et al. (1997), though they do 
not discuss in detail the magnetic field. 

This demonstrates that magnetic fields stronger than 
1 xtG exist on scales as large as the cluster virial radius 
and exert non-negligible pressure compared to the ther- 
mal one, at least around the region of the north-west 
relic. Recent observations of radio relics around other 
clusters have also revealed similarly strong magnetic fields 
of > 0.8 /xG (e.g. Hcnriksen & Mushotzky 2001; Feretti 

6 Neumann 2006; Hubert et al. 2008). That of the Abell 
3667 north-west relic is the strongest yet observed. 

It is important to note that the large magnetic field we 
have found applies only to the interior of the radio relic. 
Thus, it is possible that most clusters do not have such 
large fields at such outer radii. The relic is the brightest 
diffuse cluster radio source known, and its interior may 
not be typical of the diffuse intracluster medium outside 
the relic, or in other clusters. It may be that the mag- 
netic field and seed relativistic particles were injected into 
this region by an AGN in the past (e.g. Ensslin & Gopal- 
Krishna 2001). If the relativistic electrons were acceler- 
ated by a merger shock, the same shock may have ampli- 
fied the magnetic field, perhaps by some mechanism like 
the Weibel instability (e.g. Okabe & Hattori 2003; Fujita 
et al. 2006; Kato 2007) or the cosmic ray streaming insta- 
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bility (Lucek & Bell 2000, Bell 2004). Also, Abell 3667 
has undergone a very dramatic merger, and may not be 
typical of more relaxed clusters. 

4-3. Comparison of the HXD-PIN upper-limit of the IC 
emission with other results 

The HXD-PIN data of the NWR pointing provides an 
upper-limit on the IC emission as 5.3 x 10~ 12 erg s _1 
cm~ 2 . Since the FOV of the PIN is ~ 10 times larger 
than the region used in the XIS analysis, it provides valu- 
able limits on any non-thrcmal electrons that may be dis- 
tributed on scales larger than the north-west relic. These 
electrons, in addition, are not readily observable in the 
radio band. For example, if the relic is surrounded by 
a larger region with lower magnetic field, say 0.3 [iG, 
electrons therein with 7~5x 10 3 would be unobservable 
through their synchrotron emission since its frequency is 
less than 10 MHz. 

Using the Beppo-SAX PDS data with an equivalent ex- 
posure (113 ks) to the NWR pointing, Fusco-Femiano 
et al. (2001) derived the upper-limit hard X-ray flux of 
~ 6.4 x 10~ 12 erg s" 1 cm" 2 in the 20-80 keV band as- 
suming a r = 2.1 power-law emission. Here, they fixed 
the ICM temperature to 7 keV and included only the sta- 
tistical errors. Using the same data with similar thermal 
model, Nevalainen et al. (2004) presented a 90% confi- 
dence upper limit flux of 7.1 x 10 -12 erg s _1 cm -2 in the 
same energy band assuming a power-law with T = 2.0. 
The main difference in the two results is that the latter 
author took into account the possible systematics in the 
PDS background modeling and AGN contribution derived 
in their own manner. 

By comparing the observational results on the Crab, 
we found that the 20-80 keV flux derived from Beppo- 
SAX PDS is 21% smaller than those derived from Suzaku. 
Thus, results from Fusco-Femiano et al. (2001) and 
Nevalainen et al. (2004) can be converted into Suzaku- 
equivalent fluxes of 8.1 and 9. Ox 10~ 12 erg s _1 cm~ 2 , re- 
spectively. Therefore, our Suzaku result provides more 
stringent upper limit on the flux, compared to these two 
results. Furthermore, it is focused on a more spatially re- 
stricted region around the north-west radio relic, thanks 
to the narrow FOV of the HXD (0 - 34' FWHM) com- 
pared to that of the PDS (0-1.3° FWHM). 

4-4- Possible very hot component in the merging cluster 

The HXD-PIN results from the cluster center shows a 
"soft-shaped" hard component with T = 3.6 jig' g, which 
can be best interpreted as a very hot component with 
a temperature of 19.2tg'o keV, in addition to the — 4.7 
keV component in the 2kT picture combined with the XIS 
spectra. Using the 3kT model, the temperature of the 
very hot component is required to be higher than 20 keV. 
Since the spectra can be marginally fitted with a F ~ 1.4 
power-law component as well, the emission could be of 
non-thermal origin, or contamination from a yet unknown 
Compton-thick AGN. However, our contaminating point- 
source survey (appendix 3) and the softness of the spectral 
shape itself suggests the thermal interpretation. Since the 
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excess signal to the lkT component is significant in the 
center pointing, the very hot component can be significant 
around the cluster center. Detailed modeling of the ICM 
using all the available data, including our new XMM data 
to be obtained soon, will help clarifying the actual value. 

From the ASCA observations, this cluster is known to 
host significant temperature inhomogeniety, showing re- 
gions with a temperature as high as 10 keV together with 
the 4-6 keV dominant emission (Markevitch 1998). Later, 
Bricl, Finogucnov & Henry (2004) also showed the com- 
plex multi temperature nature of the ICM using the XMM 
data. The temperature ranges from 4 keV to 8 keV, and 
several hot regions are visible around the cluster center. 
Since these analyses are based on the projected X-ray 
color maps, the temperature of the actual hottest com- 
ponent can be higher considering the convolution along 
the line-of-sight. The Beppo-SAX PDS spectrum is con- 
sistent with the center-2fcT model taking into account all 
the systematic errors. Thus, our results on "soft-shaped" 
hard component is basically consistent with existing data. 

Such a very hot component, with a temperature as high 
as 20 keV, is suggested by many hydrodynamical simula- 
tions of merging clusters right after the time of impact 
(e.g. Takizawa 2000; Ricker & Sarazin 2001). Since the 
Abell 3667 cluster is considered to be at this stage, the 
very hot component will be the hottest heated ICM emis- 
sion. We note that a similar very hot component is sug- 
gested in the cluster RXJ1347.5-1145 (Ota et al. 2008) 
using the Suzaku wide-band data combined with Chandra 
spectra. Its high temperature and relatively large lumi- 
nosity means that the very hot component cannot be in 
pressure balance with the cooler one, otherwise, it will 
dominate the volume in the cluster. Thus, the emission 
should be temporary originating from the shock and/or 
the compression heating, and is expected to disappear 
soon by adiabatic expansion. To investigate the location 
of the emission, we will need much wider energy band cov- 
erage in the imaging detector, specifically, arcmin-scale 
imaging spectroscopy up to — 40 keV. 

5. Conclusion 

From Suzaku mapping observations of Abell 3667, 
upper-limits on the non-thermal (hard) X-ray emission 
were obtained. On the north-west radio relic itself, the 
XIS limit of 7.3xl0 -13 erg s _1 cm~ 2 extrapolated to the 
10-40 keV band provides a lower-limit magnetic field of 
> 1.6 fxG. The non-thermal energy density there is esti- 
mated to be higher than 7%, and likely near 20%. 

By combining the XIS and the XMM mapping data to 
the HXD-PIN spectra, the upper limit on the 10-40 keV 
excess hard emission is derived as 1.0-5.3x 10~ 12 erg s _1 
cm~ 2 , depending on the modeling of the ICM. Since the 
north-west relic has a magnetic field of > 1 p:G and hence 
must be almost an order of magnitude darker than this 
hard X-ray limit, the relativistic electrons must be dis- 
tributed in regions other than the relic with low magnetic 
field (of B < 0.3 /j,G) if IC hard X-ray emission near the 
flux limit exists. 
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The PIN data from the pointing on the cluster center 
shows hard X-rays in excess to the averaged ICM emission, 
which is thought to have a temperature of around 7 keV. 
Although other interpretations cannot be rejected, its soft 
spectral shape indicates the emission to be possibly from 
a very hot thermal component with a temperature of ~ 20 
keV, generated temporary in the cluster merger phase. 

In the near future, X-ray observatories with hard X-ray 
imaging optics operating up to 40 keV will be launched, 
such as the ASTRO-H (former NeXT; e.g. Takahashi 
ct al 2006), NuSTAR and Simbol-X (e.g. Fcrrando et 
al. 2006). These detectors will drastically improve the 
sensitivity to both point sources and relatively clumpy 
extended hard X-rays, such as the local non-thermal 
emission associated with the radio halos/relics and the 
very hot component. In addition, the Soft Gamma-ray 
Detector onboard ASTRO-H, (a collimated instrument) 
features the highest sensitivity at around 100 keV using 
the novel narrow-ficld-of-view Compton camera concept. 
It will open a new window to the widely distributed 
non-thermal emission, which is difficult to separate from 
the ICM thermal emission below ~ 50 keV. 
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Appendix 1. The CXB estimation for the XIS 

The Lockman Hole observation (Suzaku observation ID, 
101002010) was used as the CXB template. The NXB sub- 
tracted data from the Lockman Hole were fitted with a 
model composed of a power- law with a fixed T of 1.4, and 
thermal components (so called "Galactic components") 
with one of the temperatures fixed at 0.08 keV and an- 
other fitted to give ~ 0.3 keV. The 2-10 keV flux of the 
power-law component was derived as 5.7 ±0.1 x 10 -8 erg 
s" 1 cm -2 str -1 , which is consistent within 1 a from the 
best fit value derived using ASCA (Kushino et al. 2002). 

We also analyzed the XIS data of the outer boundary 
for Abell 3667, using the northern 1/3 region of the NWR 
pointing (hereafter the 1/3 NW- north region). It gave a ~ 
16% higher normalization for the power-law component, 
which is consistent with the expected fluctuation of CXB, 
discussed later in this section. The "Galactic component" 
parameters are similar, but the flux is brighter by a factor 
of 2-5, which is also consistent with its sky fluctuations. 
Because the 1/3 NW- north region could be contaminated 
by the outskirts of the cluster ICM emission, we modeled 
the CXB with the power-law component fixed at the value 



obtained from the Lockman Hole observation, while the 
"Galactic component" was fixed at that of the 1/3 NW- 
north region. 

The CXB fluctuations can be modeled as ctcxb /Icxb c< 
n~ - 5 S° 25 - Here, Cl e is the effective solid angle and S c 
is the upper cutoff flux. From the HEAO-1 A2 results, 
Shafer (1983) derived ct C xb//cxb = 2.8%, with n e = 15.8 
deg 2 and S c = 8 x 10" 11 erg s" 1 cm" 2 . By scaling this 
result with the XIS parameters, i.e. fi e = 0.09 deg 2 and 
S c ~lx 10" 13 erg s" 1 cm" 2 in the cluster vicinity, we 
obtain the CXB fluctuation over the XIS full FOV as 11% 
(90% confidence level). For the NWR 1/3 NW-north re- 
gion, the value is then 19%. For the north-west relic region 
with fl e = 0.04 deg 2 and 5 c ~2x 10" 14 erg s" 1 cm" 2 , it 
becomes 15.9%. Note that this scaling relation is consis- 
tent with the ASCA and Ginga CXB analysis (Ishisaki 
1997). 

Appendix 2. The CXB estimation for the HXD 

The CXB flux in the PIN is modeled as a cutoff power- 
law, described as N(E) = n x (E)~ r x exp(— E / E io \ A ) in 
photons cm" 2 s" 1 keV" 1 FOV" 1 . Here, E io \ A is an fold- 
ing energy, and n is the normalization. In energies below 
~ 80 keV, the HXD energy response for a diffuse source is 
the same to that for a point source except for its normal- 
ization. Thus, we simply scaled our CXB model to the 
solid angle of the HXD, flf XI} = 0.32 deg 2 and utilized 
the HXD nominal response. We first used the parameters 
based on the HEAO-1 results (Boldt 1987), which gave 
r = 1.29, E iold = 40.0 and n = 8.36 x 10" 4 . We then com- 
pared this model to the XIS observation of Lockman Hole, 
on which our XIS CXB modeling is based on. It was found 
that the model predicts a flux 8.6% lower in the XIS 3-8 
keV band data. Thus, we modified the PIN CXB model 
normalization to match the difference, as n = 8.69 x 10" 6 . 
Note that although there is an apparent difference in pho- 
ton index between the CXB model of the XIS and PIN, 
the cutoff function in the later makes two models perfectly 
matched in the energy band at 3-8 keV. Through the same 
approach applied to the XIS data, the CXB fluctuation in 
the HXD FOV is expected to be 18%, with f^ XD = 0.32 
deg 2 and a conservative upper cutoff flux of S c ~ 8 x 10" 12 
erg s _1 cm" 2 in the 10-40 keV band. 

Appendix 3. Estimation of point source contam- 
ination in the FOV of the PIN 

Since the PIN has a wider FOV (68' x 68' bottom to bot- 
tom) than that of the XIS (18' x 18') and also lacks imag- 
ing capability, signal contamination from point sources 
within the FOV should be carefully estimated. Figure 13 
shows the PSPC image of Abell 3667 with the FOVs of the 
PIN superposed. The PSPC data is processed through the 
method described by Snowden et al. (1994). There exists 
seven X-ray point sources in the PIN FOVs. We identified 
these sources using the SIMBAD astronomical database. 

Source 4 and Source 6 are in the FOV of XIS, and thus 
their flux is estimated by the XIS data. Although the 
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Fig. 13. The background subtracted ROSAT PSPC image 
(50' radius) of Abell 3667. Gray-scale contour is logarithmic 
scale. The FOVs (68' X 68' bottom to bottom) of PIN are 
shown in dashed thin rectangles, those of the XIS are shown 
in solid rectangles, and X-ray point sources in the FOVs are 
shown in circles. 

other sources are outside the FOV of XIS, archival data of 
ASCA GIS are available for the sources 1, 2, 3, and 5. We 
therefore extracted the GIS spectra of source 1 thorough 
source 6 and fitted them with a power low modified by 
galactic absorption. Blank sky observations of the GIS 
were used for the background of source 1,2 and 3, while 
neighboring field is used for source 5. 

Using the best-fit parameters, we simulated the PIN 
spectra of the six sources for the three pointings, consid- 
ering the PIN transmission to each of them on the three 
pointings. The total contribution of the six point sources 
thus estimated is more than an order of magnitude lower 
than the detected signal, and hence negligible. Source 7, 
not included in this estimation, is located at the boundary 
of the PIN FOV of the center pointing. It is also a weak 
hard X-ray source, since it is a distant cluster of galaxies. 
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